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Peculiarities of obtaining parameters for broad multi-channel resonances from data are discussed 
analyzing the experimental data on processes ittt — > tytt,KK in the I G J PC = + ++ channel 
in a model-independent approach based on analyticity and unitarity and using an uniformization 
procedure. We show that it is possible to obtain a good description of the tx-k scattering data from the 
threshold to 1.89 GeV with parameters of resonances cited in the PDG tables as preferred. However, 
in this case, first, representation of the tttv background is unsatisfactory; second, the data on the 
coupled process rnr — > KK are not well described even qualitatively above 1.15 GeV when using the 
resonance parameters from the only Tin scattering analysis. The combined analysis of these coupled 
processes is needed, which is carried out satisfactorily. Then both above-indicated flaws, related 
to the analysis of solely the 7T7r-scattering, are cured. The most remarkable change of parameters 
with respect to the values of only tt-k scattering analysis appears for the mass of the /o(600) which 
is now in some accordance with the Weinberg prediction on the basis of mended symmetry and 
with an analysis using the large- N c consistency conditions between the unitarization and resonance 
saturation. The obtained 7T7r-scattering length a{j in case when we restrict to the analysis of the rnv 
scattering or consider so-called A-solution (with a lower mass and width of /o(600) meson) agrees 
well with prediction of chiral perturbation theory (ChPT) and with data extracted at CERN by the 
NA48/2 Collaboration from the analysis of the K e 4 decay and by the DIRAC Collaboration from 
the measurement of the -k + tt~ lifetime. 
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I. INTRODUCTION 



The study of scalar mesons is very important for understanding the QCD vacuum. However, despite of the big 
effort devoted to studying various aspects of the problem [l[ (for recent reviews see, e.g. 043) a description of this 
mesonic sector is far from being complete. Parameters of the scalar mesons, their nature and status of some of them 
are still not well settled [l[ . For example, applying our model- independent method in the three-channel analyses of 
processes 7T7t — > nir,KK, and r\r\ or 7777' 6, 7] we have obtained parameters of the /o(600) and /o(1500) which differ 
considerably from results of analyses which utilize other methods (mainly those based on the dispersion relations 
and Breit-Wigner approaches). This difference is very interesting and its reasons should be understood because our 
method of analysis is based only on the demand for analyticity and unitarity of the amplitude using a uniformization 
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procedure. The construction of the amplitude is practically free from any dynamical (model) assumptions utilizing 
only the mathematical fact that a local behaviour of analytic functions determined on the Riemann surface is governed 
by the nearest singularities on all corresponding sheets. I.e., the obtained parameters of resonances can be considered 
as free from theoretical prejudice. 

To better understand reasons for the difference in results, in this paper we have performed first the two-channel 
model-independent analysis only of the tttt scattering data. The thing is that in our previous three-channel analysis 
with the uniformizing variable [6, 7] we were enforced to construct a four-sheeted model of the initial eight-sheeted 
Riemann surface. This we have achieved by neglecting the lowest 7T7r-threshold branch-point which means that we 
have considered the nearest to the physical region semi-sheets of the initial Riemann surface. This is in the line with 
our approach of a consistent account of the nearest singularities on all relevant sheets. The two-channel analysis 
utilizes the full Riemann surface and is, therefore, free of these approximations. 

In the two-channel analysis only of the 7T7r scattering data we have obtained a good description from the threshold 
to 1.89 GeV with values of parameters of the /o resonances, which are well consistent with the ones cited in the PDG 
tables [H as preferred. However, it was turned out that the cross-section of the coupled process tttt — » KK is not well 
described even qualitatively above 1.15 GeV when using for the relevant resonances the same values of parameters. 
The combined analysis of these coupled processes is needed, which is carried out also satisfactorily. 

In the presented combined two-channel analysis of data on processes mr mr and tttt —> KK, we have checked 
whether the results of our previous three-channel analysis 0, 0] are also obtained in the two-channel case and therefore 
shown whether the above-indicated assumptions are justified. It is also interesting to determine and discuss the 
scattering length ajj and slope bg, related to the effective range of interaction, in the separate analyses using the 
alternative data on tttt scattering from Refs. [|[ and ||, adding the very precise NA48/2-Collaboration 7T7r-data [l(j. 

In Sect. II a basic formalism for our two-channel model- independent method is shown and rules for the calculation 
of resonance parameters discussed. Results of analyses are presented and discussed in Sect. Ill, first only for the tttt 
data to clarify a consistency of our results with values of parameters from the PDG tables and then for the combined 
two-channel analysis, using two sets of data for the 7T7r scattering [H Q and all accessible data 11 1 for tttt — > KK, to 



verify a plausibility of our assumptions in the three-channel calculations. Conclusions are provided in Sect. IV. 



II. THE COUPLED-CHANNEL FORMALISM IN UNIFORMIZING VARIABLE METHOD 

Our "model-independent" method which essentially utilizes an uniformizing variable can be used without any 
further assumptions only for the two-channel case. In the three-channel case, some assumptions about the Rieman 
surface have to be made @, 0). In this work we consider the two-channel case. 

The two-channel S'-matrix is determined on the four-sheeted Riemann surface. The matrix elements iSy, where 
i,j = 1(tttt),2(KK) denote channels, have the right-hand cuts along the real axis of the s complex plane (s is the 
invariant total energy squared), starting with the channel thresholds Si (i = 1,2), and the left-hand cuts related to 
the crossed channels. The Ricmann-surface sheets are numbered according to the signs of analytic continuations of 

the roots yjs — Sj (i — 1,2) as follows: signs(lm-y/s — s^Imys — 52)= ++1 — h >H — correspond to sheets I, II, 

III, IV. 

The resonance representations on the Riemann surface are obtained from the formulas which express an analytic 
continuation of the S'-matrix elements to unphysical sheets in terms of the matrix elements on sheet I (physical) 
having only resonance zeros (beyond the real axis), at least, around the physical region: 
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Then, starting from the resonance zeros on sheet I, one can obtain an arrangement of poles and zeros of a resonance 
on the whole Riemann surface. 

In the one-channel consideration of the elastic scattering 1 — > 1 the main model-independent contribution of a 
resonance is given by a pair of conjugate poles on sheet II and by a pair of conjugate zeros on sheet I at the same 
conjugate points of complex energy in S11. In the two-channel consideration of the processes 1—^1 and 1 — > 2, a 
resonance is represented by a pair of conjugate poles on sheet II and by a pair of conjugate zeros on sheet I in Sn and 
also (as it is seen from eq.([2])) by a pair of conjugate poles on sheet III and by a pair of conjugate zeros on sheet IV 
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at the same conjugate points of complex energy if the coupling of channels is absent (512 = 0). If a resonance decays 
into both channels and/or takes part in processes of exchange in the crossing channels, the coupling of channels arises 
(S12 7^ 0). Then positions of the conjugate poles on sheet III (and of corresponding zeros on sheet IV) turn out to 
be shifted with respect to the positions of the zeros on sheet I. Thus we obtain a cluster of poles and zeros (the pole 
cluster of type (a)) which gives the main mo del- independent contribution of the corresponding two-channel resonance. 
Obviously, depending on nature of a resonance there are two more pole clusters: when the pair of conjugate zeros on 
sheet I, corresponding to the resonance, is present only in S22 - the pole cluster of type (b) - and when in each of 
Sn and £22 ~ °f type (c). For the resonances of type (b), the pair of complex conjugate poles on sheet III is shifted 
relative to the pair of poles on sheet IV. For the states of type (c) , one must consider the corresponding two pairs of 
conjugate poles on sheet III. 

It is clear that for calculating the resonance parameters (masses, total widths and coupling constants with channels) 
one must use the poles on those sheets where they are not shifted due to the coupling of channels as they respects 
the zero position on the physical (I) sheet. For resonances of types (a) and (b) these poles are on sheets II and IV, 
respectively. For resonance of type (c) the poles can be used on both these sheets. 

Analogous consideration can be carried out in the three-channel case Seven types of resonances arise there. 

Formulas of analytic continuations of the S-matrix elements to unphysical sheets in the general case of N channels can 
be found in [!§]• In these cases one can see that only on the sheets with the numbers 2 l (i = 1, • • • , N is the number 
of channel), i.e. II, IV, VIII,. . ., the analytic continuations of the S'-matrix elements have the form oc where 
S}i is the S'-matrix elements on the physical (I) sheet. This means that the pole positions of resonances only on these 
sheets are at the same points of the s-plane, as the resonance zeros on the physical sheet, and are not shifted due 
to the coupling of channels. Therefore, the resonance parameters should be calculated from the pole positions only on 
these sheets. It seems that neglecting this fact serious misunderstandings can arise. This concerns analyses which do 
not consider the structure of the Riemann surface of the S-matrix and especially the analyses of the decay processes 
in which, as a rule, the multi-channel nature of resonances is not taken into account. 

In the combined analysis of coupled processes, it is convenient to use the Le Couteur-Newton relations [l4j]. They ex- 
press the S-matrix elements of all coupled processes in terms of the Jost matrix determinant d(y/s — si, • • • , y/s — sn) 
that is a real analytic function with the only square-root branch-points at y/s — s a = 0. 

A necessary and sufficient condition for existence of the multi-channel resonance is its representation by one of the 
types of pole clusters. In order to use the representation of resonances by the pole clusters which is very important 
for the wide multi-channel states, we take advantage of the fact that the amplitude is a one-valued function on the 
Riemann surface. To this end, a uniformizing variable is applied, which maps the Riemann surface onto a plane. 

In this two-channel analysis of processes tttt tttt, KK we applied the uniformizing variable which takes into 
account, in addition to the tttt- and ifif-threshold branch-points, the left-hand branch-point at s = 0, related to the 
tttt crossed channels: 

v _ rn K y/s- 4ml + m 7r y/s - Am 2 K 
y/s(m 2 K - m\) 

It maps the four-sheeted Riemann surface with two unitary cuts and the left-hand cut onto the u-plane [l5j | divided 
into two parts by a unit circle centered at the origin. In Figure UH an uniformization u-plane for the two-channel- tttt- 
scattering amplitude is shown with the representation of resonances of types (a), (b) and (c): the Roman numerals 
(I,. . . , IV) denote the images of the corresponding sheets; the thick line represents the physical region; the points 
i, 1 and b = yj (juk + rn^)/ (tok — m^) correspond to the tttt, KK thresholds and s = 00, respectively; the shaded 
intervals (-co, —b], [— £> -1 , b^ 1 ], [6, 00) are the images of the corresponding edges of the left-hand cut. The depicted 
positions of poles (+) and of zeros (o) give the resonance representations of the type (a), (b) and (c) in Sn. The 
resonance poles are symmetric to the corresponding zeros with respect to the unit circle that guarantees the elastic 
unitarity of tttt scattering up to the KK threshold. The whole picture of poles and zeros is symmetric with respect 
to the imaginary axis guaranteing a property of the real analyticity of the S-matrix. 

On the w-plane, Sn(y) has no cuts; Si2(v) and S22(y) do have the cuts which arise from the left-hand cut on the 
s-plane, starting at s = 4{m 2 K — ml), which is further approximated by a pole 

d L = v- i (l - (p- i\J\ - F» 4 (l + [p + i\J\ - p 2 )v) 4 (from analysis p = 0.903 ± 0.0004). (3) 

The fourth power is stipulated by the following model- independent arguments. First, a pole on the real s-axis on the 
physical sheet in S22 is accompanied by a pole on sheet II at the same s-value (see eq.Q). On the w-plane this implies 
the pole of second order (and also zero of the same order, symmetric to the pole with respect to the real axis). Second, 
for the s-channel process tttt — > KK, the crossing u- and i-channels are the tt — K and W — K scattering (exchanges 
in these channels give contributions on the left-hand cut). This results in the additional doubling of the multiplicity 
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FIG. 1: Representation of resonances of type (a), (b) and (c) on the uniformization u-plane in Sn. 



of the indicated pole on the i>-plane. Therefore, the contribution of the left-hand branch-point at ,s = A(m 2 K — m 2 ) is 
approximated as the fourth-power pole on the real s-axis on the physical sheet in the sub-ifA'-thrcshold region. 
On the u-plane, the Le Couteur-Newton relations are [H], [H| 

_ d(-v^) d(v-i) d(-v) 

11 ~ ~ ' 22 6llA22 ~ 5l2 -^T (4) 

The main model-independent contribution of resonances, given by the pole clusters, is factorized in the ^-matrix 
elements from the background. The possible remaining corrected and model-dependent contributions of resonances 
are supposed to be included in the background. Therefore the d(w)-function, which does not possess already branch 
points, is taken as 

d = d res d L dbg. (5) 
The function d res (v) represents the contribution of resonances, described by one of three types of the pole-zero clusters, 



i.e., 

M 



-M 



H(l-V*v)(l+VnV), (6) 



where M is the number of pairs of the conjugate zeros. 
The background part is 



d bg = exp[-i ^ Sn (a n + ifin)] (7) 



71=1 
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with 



Oral 



(s - s a ) + a n 



Sin = Ki + h nv -9(s - s v ) + b na -6(s - s a ) + b nv -9(s - s v ) 

Sn &<t °v 

where and s a are the r\r\ and 00 thresholds, respectively (the latter should be determined in the analysis), s v is a 
combined threshold of many opened channels in the vicinity of 1.5 GeV (e.g., 7777', pp, ujlo) and it is determined in 
the analysis: s CT = 1.6558 GeV 2 , s v = 2.1293 GeV 2 . 

The data used in the analysis are the results of phase analyses which are given for phase shifts of the amplitudes 
8 a p and for modules of the S'-matrix elements 77 Q( 3 = \S a p\ {ot.,fi = 1, 2) 

S aa = v aa e 2lS - , S af} = tiape*-* . (8) 

The two-channel unitarity provides the relations 

Vn = V12 = (1 - ^ii 2 ) 1/2 ; 012 = Su + 5 2 2- (9) 



III. ANALYSIS OF THE DATA ON ISOSCALAR S-WAVE PROCESSES nn -> nn, KK 

We analyzed data on processes tttt —> tttt,KK. For the 7T7r scattering, we took alternative data - set I: for 
0.575 GeV < y/s < 1.89 GeV from Ref. @| and for ^ < 1 GeV from Refs. jlfl, MM, M, H3 ; set II: for 0.61 GeV < 
y/s < 1.59 GeV from Refs. @ and [T^ (solution G) and for y/s < 1 GeV as in set I. For -kit —> KK, we used practically 
all accessible data (Tlj . 

Initially analyzing only the tttt scattering from set I, we demonstrated that it is possible to achieve an excel- 
lent description of the data for the phase shift 6u and modulus of the S'-matrix element (the total x 2 /NDF = 
168.433/(189— 33) s» 1.08) with the parameters of resonances (TableU)) which largely coincide with the values cited 
as estimation of the PDG though a negative phase-shift in the background on the 7T7t threshold arises. E.g., 
for the Jo (600) the found pole on sheet II coincides practically with the one at around 450 — i275 MeV which was 
found in the recent dispersive 7T7r-scattering data analyses [26l |27| . More remarkable differences are only for the mass 
of / (980) (1001 MeV against 980±10 MeV of PDG) and the width of / (150Q) (336 MeV against 109±7 MeV of 
PDG). However, the mass of /o(980) slightly above 1 GeV was also obtained in many other works which analyzed 
the 7T7r scattering (e.g. [H|)- As to the latter, we think that the observed wide resonance /o(1500) is in reality a 
superposition of two states, wide and narrow. The narrow state is observed in processes considered in works cited by 
the PDG. To test this interpretation of the /o(1500), we analyzed the 7T7r-scattering data also assuming both wide 
and narrow / (1500). The gained description is excellent: the total x 2 /NDF = 171.715/(189 - 29) w 1.07. The 
obtained parameters of resonances are shown in Table |TJ Now the parameters of the narrow /o(1500) are consistent 
with those in the PDG tables. In the presented analyses, the /o(600) and /o(980) are described by the clusters of 
type (a); / (1370), /o(1500) and /o(1710), type (b); /q(1500), type (c). The received background parameters in the 
analysis with the narrow / (1500) are: an = -0.0895 ± 0.0030, a ln = 0.04 ± 0.03, a la = 0.0 ± 0.8, a lv = 0.0 ± 0.7, 
611 = 0.0 ± 0.007, bi v = 0.0 ± 0.01, b la = 0.0 ± 0.02, b lv = 0.054 ± 0.036. 

In Figure [2] results of the fitting only to the 7T7r-scattering data are shown (upper row); in the lower row there are 
given energy behaviours of the phase shift and module of the irir —} KK matrix element which are calculated using 
the resonance parameters from the analysis of only 7T7r-scattering: the dotted and short-dashed lines correspond to 
the analysis without and with the narrow /o(1500), respectively. The gained description of data is very good. The irir 
scattering length a®, obtained in these analyses of the data from set I, is 0.222 ± 0.008 m~\ , which is also in the very 
good agreement with the experimental results and with the results of the ChPT calculations (see Tab. V). However, 
let us also emphasize two important flaws: 

• First, the negative phase-shift in the background beginning from the irir threshold (an = —0.0895) seems to 
be necessary for a sucessful description of the data. This should not be the case because, in the uniformizing 
variable, we have taken into account the left-hand branch-point at s = which gives a main contribution to 
the 7T7T background below the KK threshold. Other possible contributions of the left-hand cut related with 



exchanges by the nearest mesons - the p-meson and the /o(600) - practically obliterate each other [15j because 
vector and scalar particles contribute with the opposite signs due to the gauge invariance. 
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TABLE I: Pole clusters for resonances on the complex energy plane y/s in the analysis of only 7T7r-scattering without and with 
the narrow /o(1500). Pole energies ^s7 = E r — ir r /2 in MeV are shown. 



Without the narrow /o(1500) 


Sheet 


II 


III 


IV 


/o(600) 


r r /2 


A A ^7 1 f A 

447.5±5.4 
267.0±6.2 


492.6±31.5 
307.8±15.6 




/o(980) 


r r /2 


1001. 1±3. 4 
20.3±2.4 


978. 5± 9.5 
38.5± 6.3 




/o(1370) 


E r 




1382.6±38.5 
179.5±39.2 


1301.2±38.2 
243.0±52.8 


/o(1500) 


E r 

Tr/2 


1512.0±10.6 
191.0±10.5 


1499.0±108.3 1509.1±112.6 
310.2± 71.8 241.0±68.9 


1505.9±30.0 
168.0±32.8 


/o(1710) 


E r 

r r /2 




1700.3±30.5 
58.8±29.5 


1720.1±30.5 
64.9±34.3 


With the narrow /o(1500) 


Sheet 


II 


III 


IV 


/o(600) 


E r 
Tr/2 


447.5±5.9 
267.0±6.5 


492.7±36.0 
307.8±16.5 




/o(980) 


E r 

Tr/2 


1001.1±3.7 
20.3±2.6 


979.1±12.0 
38.5±7.1 




/o(1370) 


E r 
Tr/2 




1375.8±51.5 
179.5±36.5 


1301.1±47.9 
224.0±49.3 


/o(1500) 


E r 

Tr/2 




1498.8±39.3 
51.8±43.3 


1503.7±45.1 
56.5±39.4 


/o(1500) 


Er 
Tr/2 


1511.4±11.2 
200.5±11.0 


1499.8±104.3 1509.1±119.4 
310.5±58.8 241.0±63.8 


1505.9±38.5 
168.1±40.6 


/o(1710) 


E r 

r r /2 




1700.3±31.2 
58.6±26.4 


1720.1±32.2 
64.9±30.1 



• Second, the description of the data on reaction tttt — > KK 1 using the same parameters of resonances as in 
the 7T7r channel, is satisfactory only for the phase shift 4>i2 which is due to the approximation of the left-hand 
branch-point at s = A{m 2 K — m 2 ) in S±2 and S22 by the fourth-power pole. The module of the S-matrix element 
7712 is described well only from the KK threshold upto about 1.15 GeV as it should be due to the two-channel 
unitarity (see eqs. ©). Above this energy the description fails even qualitatively (see Fig. [2]). 

From this we conclude that: If the data are consistent, for obtaining correct parameters of wide resonances the 
combined analysis of data on coupled processes is needed. Further the combined analyses of data (sets I and II) on 
processes tttt — > tttt,KK are performed supposing that in the 1500-MeV region there are two resonances. 

In the analysis of set I, the resonances are described by pole clusters of the same types as in the analysis only 
of 7T7r-scattering. Satisfactory combined description of two analyzed processes is obtained with the total x 2 /NDF = 
391.299/(312 — 40) ss 1.44. One sees that the data for the tttt scattering below 1 GeV admit two solutions for the phase 
shift: A and B which mutually differ mainly in the pole position on sheet II for the /o(600). The \ 2 shown above is 
for the solution B. The A-solution gives a slightly worse result: the total x 2 / NDF = 416.887/(312 - 40) « 1.53. 

In the analysis of set II, the resonances are described by the pole clusters of the same types as in previous case 
except for the /o(1500) which is represented now by the cluster of type (a) to provide more rapid growth of the 7T7t 
phase shift above 1.45 GeV than in the first case. 

Also in this case, a satisfactory description is obtained with the total x 2 /NDF = 418.268/(306— 41) w 1.58 for the 
A-solution and xV NDF = 375.160/(306-41) 1.42 for the B-solution. 

In Figure |3] results of the fitting to the experimental data from sets II are shown. The solid lines correspond to the 
A-solution and the dashed ones to the B-solution. 

In Table [TT] we show the pole clusters for resonances on the complex energy plane yfs for the A- and B-solutions 



7 




FIG. 2: The S-wave phase shifts and modules of the 7T7r-scattering matrix element (upper row) and of the nn — > KK matrix 
element in analyses of set I. The dotted and short-dashed lines correspond to the analysis of only 7T7r-scattering without and 
with the narrow /o(1500), respectively. The long-dashed and solid lines correspond to solutions A and B of the combined 
analysis of tttv — > tttt,KK, respectively. The data are from Refs. [1. HoL UlL IT^ - Ha . l20l. |2~D] . 



in the combined analyses of data (set I) on processes tttt — > 7T7t, KK. The obtained background parameters in 
the analysis of set I for the A-solution are: an = 0.0 ± 0.003, a lr) = -0.1004 ± 0.0301, a 1<7 = 0.2148 ± 0.0822, 
axv = 0.0 ± 0.07, b n = b Vq = b Xa = 0, b lv = 0.012 ± 0.0287, a 21 = -0.919 ± 0.107, a 2r) = -1.399 ± 0.348, 
a 2a = 0.0 ± 0.7, a 2v = -11.45 ± 0.75, b 2 i = 0.0747 ± 0.0503, b 2ri = b 2<J = 0, b 2v = 4.83 ± 1.94; 

for B-solution: a n = 0.0±0.003, a lv = -0.0913±0.0327, a la = 0.1707±0.0899, a lv = 0.0±0.07, b n = b lv = b la = 0, 
b lv = 0.006±0.029, a 21 = -1.338±0.111, a 2v = -I.119i0.376, a 2a = 0.0i0.8, a 2v = -12.13i0.77, b 21 = 0.018i0.050, 
b 2v = b 2a = 0, b 2v = 4.48 i 1.98. 

It is apparent that in the combined analysis of data on coupled processes both above-indicated important flaws, 
which related to the analysis of only 7T7r-scattering, are cured. Now the 7T7t background below the KK threshold is 
absent (an = 0.0) because its contribution is practically completely accounted for by the left-hand branch-point at 
s = which is included explicitly in the uniformizing variable @. An arising pseudo-background at the 7777-threshold 
(ax?] < 0) is also clear: this is a direct indication to consider explicitly the 7777-threshold branch-point. This was already 
done in our previous work Q- 

Generally, wide multi-channel states are most adequately represented by pole clusters, i.e., by the poles on all 
the corresponding sheets, because the pole clusters give a main effect of resonances. The pole positions are rather 
stable characteristics for various models, whereas masses and widths are very model-dependent for wide resonances 
(see a discussion in Ref. [IB])- Earlier one noted that the wide resonance parameters are largely controlled by the 
nonresonant background (see, e.g. [23I]). In part this problem is removed due to allowing for the left-hand branch- 
point at s = in the uniformizing variable. There remains only a considerable dependence of resonance masses and 
widths on the used model. E.g., if for the resonance part of the amplitude one use the form 

Tres = VsTel 

m res ~ S~ ly/sYtot 

then masses and total widths (Tables IIIII and IIV[) can be calculated using formulas 



= \/ E " -■- ( y) and T tot = r r . (11) 
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FIG. 3: The S-wave phase shifts and modules of the 7T7r-scattering matrix element (upper row) and of the irn — > KK matrix 
element in the combined analyses of the data on these two processes from set II. The solid lines correspond to the A-solution 
and the dashed ones to the B-solution. The data are from Refs. [jj-Ull [lr3 - [2ll | . 



where pole parameters ^fs^ = E r — ir,./2must correspond to the pole positions on sheets II and IV for the resonances 
of type (a), (c) and (b), respectively. 

Masses and total widths of the resonances, obtained in the analyses of both sets of data (Tables Hill and ITV |) . are 
reasonably close each other taking into account their errors. 

In Table W\ we compare our results for the 7T7t scattering length Oq, obtained in the analyses of the data of sets I 
and II, with results of some other theoretical and experimental works. 

Let us emphasize that in the case when we limit ourselves to the analysis of the 7T7t scattering and in the case of 
the scenario A (lower mass and width of fo (600) meson) we reproduce with a high accuracy the results of the chiral 
perturbation theory (ChPT) [2^, [2(| including constraints imposed by the Roy's equations. On the other side, the 
solutions of the scenario B (with a heavier mass and width of /o(600) meson) is similar to the the predictions of 
the chiral approaches based on the linear realization of chiral symmetry (models of the Nambu-Jona-Lasinio (NJL) 
type [H,[l3|) Taking into account very precise experiments at CERN performed by the NA48/2 Collaboration [24 



and the DIRAC Collaboration 3l|, which confirmed the prediction of the ChPT [25j, |26j we prefer the A-solution. In 



particular, the NA48/2 Collaboration [24[ extracted the S'-wave 7T7t scattering lengths 

a° = (0.2220 ± 0.0128 stat ± 0.0050 syst ± 0.0037 th ) m^ 1 

a a 2 = (-0.0432 ± 0.0086 sta t ± 0.0034 syst ± 0.0028 th ) m~\ (12) 
from the analysis of the K e 4 decay K ± — > TT + TT~e ± v. The DIRAC Collaboration extracted the quantity 



n ocqo+0.0080 
u.zooo_ 007g 



+0.0078 




-0.0073 




stat 


syst/ 



m-l (13) 



from the measurement of the 7r + 7r atom lifetime r in the ground state using the model- independent formula derived 
in Refs. [HJ at next-to-leading order (NLO) in isospin breaking: 



r 



- 1 ^(a°- a °) 2 (l + <5), (14) 
where the quantity 8 = (5.8 ± 1.2) x 10 -2 encodes the NLO isospin-breaking correction. 
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TABLE II: Pole clusters for resonances on the ^/i-plane in the combined analyses of data (set I) on processes tv-k — > -kit, KK. 
s /s r ~='E r — iTr/2 in MeV are shown. 



A-solution 


Sheet 


II 


III 


IV 


/o(600) 


rv/2 


517.0±7.8 
393.9±6.0 


458.5±14.7 
205.9±4.7 




/o(980) 


r r /2 


1004.6±3.9 
25.0±2.3 


995.5±10.1 
96.9±2.7 




/o(1370) 


E r 

Tr/2 




1351.5±32.5 
369.0±45.7 


1342.9±12.2 
221.6±30.7 


/o(1500) 


E r 

Tr/2 




1498.7±5.8 
56.7±5.6 


1501.1±6.4 
56.6±6.0 


/o(1500) 


E r 
T r /2 


1532.2±12.4 
323.2±21.0 


1489.1±16.2 1515.9±29.2 
217.9±10.2 388.4±22.6 


1519.3±18.7 
339.5±42.2 


/o(1710) 


E r 

Tr/2 




1701.9±31.8 
77.8±18.0 


1717.0±34.9 
72.9±16.2 


B-solution 


Sheet 


II 


III 


IV 


/o(600) 


E r 

Tr/2 


550.6±9.0 
502.1±7.2 


664.5±12.1 
188.2±2.6 




/o(980) 


E r 

Tr/2 


1003.2±3.0 
28.9±2.0 


995.4±7.3 
96.7±2.7 




/o(1370) 


E r 
Tr/2 




1353.8±27.9 
367.4±37.4 


1336.7±14.1 
251.9±27.5 


/o(1500) 


E r 

Tr/2 




1499.5±6.0 
56.5±6.1 


1500.3±6.3 
57.0±6.4 


/o(1500) 


E r 
Tr/2 


1528.4±12.5 
328.0±20.2 


1491.3±15.8 1510.8±29.1 
217.9±8.0 388.3±16.3 


1515.6±17.0 
340.3±34.9 


/o(1710) 


E r 

Tr/2 




1703.1±31.5 
81.7±19.9 


1722.0±35.7 
92.3±20.3 



In Table IVTl we show results for the slope parameter 6q defined in following expansion around the threshold 27 1 



^ sin2<Ju(a) = a° k + b°k 3 +(D(k 5 ), (15) 



4to„ 



where k = s /4 — is the pion cm. momentum. Our results agree well with the other results only in the case 

of the one-channel analysis (only 7T7t scattering). In the combined analysis the obtained values are by 20-30% smaller 
than the results consistent with the dispersion relations. The results do not differ so much for the A and B solutions 
as the scattering length. 

For convenience of usage of our results, we show the positions of zeros on the u-plane, which correspond to the 
resonances, obtained for the solutions A: for set I - 
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TABLE III: Masses and total widths of the resonances, obtained in the analysis of set I. 





A-solution 


B-solution 


State 


m res [MeVJ 


1 i 0t [MeVJ 


m res [MeVJ 


1 t0 i[MeVJ 


/o(600) 


650.0±9.8 


787.8±12.0 


745.2±11.5 


1004.2±14.4 


/o(980) 


1004.9±4.0 


50.0±4.6 


1003.6±3.1 


57.8±4.0 


/o(1370) 


1361.1±17.0 


443.2±61.4 


1360.2±18.9 


503.8 ±55.0 


/o(1500) 


1502. 2±6. 6 


113.2±12.0 


1501.4±6.5 


114.0±12.8 


/o(1500) 


1565.9±16.5 


646.4±42.0 


1563.2±16.8 


656.0±40.4 


/o(mo) 


1718.6±35.6 


145.8±32.4 


1724.5±36.7 


184.6±40.6 



TABLE IV: Masses and total widths of the resonances, obtained in the analysis of set II. 





A-solution 


B-solution 


State 


m res [MeV] 


r tot [MeV] 


m res [MeV] 


r iot [MeV] 


/o(600) 


652.3±11.8 


828.6±14.0 


743.4±14.1 


1001.8±17.4 


/o(980) 


1004.4±1.8 


54.4±3.6 


1004.0±1.9 


56.4±3.8 


/„(1370) 


1367.8±13.5 


444.8±53.8 


1350.1±38.0 


491.4±121.6 


/o(1500) 


1498.7±5.3 


113.6±9.8 


1499.4±5.3 


115.2±9.8 


/o(1500) 


1563.8±30.1 


656.8±75.6 


1564.1±29.7 


658.8±68.4 


/o(1710) 


1721.1±73.2 


223.8±152.4 


1721.3±47.1 


223.8±105.0 



TABLE V: The tttt scattering length • 



a° [m-l] 


Remarks 


References 


0.222 ± 0.008 


Analysis only of tttt scattering 


This paper 


0.230 ±0.004 
0.282 ±0.003 
0.226 ± 0.004 
0.275 ±0.004 


A-solution, set I 
B-solution, set I 
A-solution, set II 
B-solution, set II 


This paper: combined analysis of 
processes mr — > irn, KK 


0.26 ±0.05 


Analysis of the K — > nnev 
using Roy's equation 


L.Rosselet et al. [20] 


0.24 ±0.09 


Analysis of ir~p — > -K + -R~n 
using the effective range formula 


A.A.Bel'kov et al. [21] 


0.2220 ± 0.0128 stat ± 0.0050 syst ± 0.0037 th 


Experiment on K e 4 decay 


J.R. Batley et al. [24] 


0.220 ±0.005 


ChPT + Roy's equations 


G. Colangelo et al. [25, 26] 


0.220 ±0.008 


Dispersion relations and K e 4 data 


R.Garci'a-Martm et al. [27] 


0.26 


NJL model 


M.K.Volkov [28] 


0.28 


Variant of the NJL model 


A.N. Ivanov and N. Troitskaya [29] 
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TABLE VI: The -kit scattering slope parameter 



b°o [m~l] 


Remarks 


References 


0.295 ±0.021 


Analysis only of irn scattering 


This paper 


0.210 ±0.010 
0.201 ±0.007 
0.209 ±0.011 
0.208 ±0.011 


A-solution, set I 
B-solution, set I 
A-solution, set II 
B-solution, set II 


This paper: combined analysis 
of processes tttt — > Tin, KK 


0.278 ± 0.005 


Analysis using Roy-like equations 
and forward dispersion relations 


Garcia-Martin et al. [2JJ 


0.290 ± 0.006 


Analysis using Roy's equations 
and forward dispersion relations 


Kamiriski et al. [35J 



for / (600) 


V\ = 


1.3600 ±0.0080 4 


-(0.3797 ±0.0076)?', 


V2 


= 0.6660 ±0 


0130- 


(0.3254 ±0.0079)?, 


for Jo (980) 


v[ = 


1.0657 ± 0.0031 4 


- (0.0346 ± 0.0023)?', 


v' 2 


= 0.9006 ±0 


0019- 


(0.0725 ±0.0034)?;, 


for / (1370) 


V3 = 


-1.2331 ±0.0010 


- (0.0419 ±0.0023)?, 




= 0.7966 ±0 


0011- 


(0.0374 ± 0.0020)?', 


for / (1500) 


v' 3 = 


-1.2467 ±0.0007 


- (0.0077 ±0.0006)?, 


< 


= 0.8023 ±0 


0004- 


(0.0050 ±0.0003)?, 


for #(1500) 


V5 = 


1.2641 ±0.0014 4 


- (0.0362 ±0.0016)?, 




= -1.2643 ± 


0.0032 


- (0.0386 ± 0.0045)? 




Vj = 


0.7981 ±0.0015- 


- (0.0184 ±0.0012)?;, 


V8 


= 0.7877 ±0 


0033- 


(0.0265 ±0.0031)?;, 


for / (1710) 


v" - 

v 3 — 


-1.2701 ±0.0030 


- (0.0062 ±0.0011)?;, 


v'l 


= 0.7881 ±0 


.0018- 


(0.0042 ± 0.0007)?; 



for set II - 



for Jo (600) 


vi = 


1.3815 ± 0.0102 ± 


(0.3713 ± 0.0096)?;, 


''2 


= 0.6360 ±0 


0159- 


(0.3286 ±0.0089)?\ 


for Jo (980) 


v[ = 


1.0664 ± 0.0022 ± 


(0.0373 ± 0.0021)?;, 


v' 2 


= 0.8998 ±0 


0017- 


(0.0742 ± 0.0025)?, 


for / (1370) 


V3 = 


-1. 2342 ±0.0013 - 


- (0.0413 ± 0.0027)?;, 




= 0.7970 ±0 


0012- 


(0.0393 ± 0.0025)?', 


for / (1500) 


v' 3 = 


1.2463 ±0.0006 - 


(0.0078 ± 0.0006)?, 




= 0.8029 ±0 


0003- 


(0.0054 ±0.0003)?, 


for #(1500) 


V5 = 


1.2642 ± 0.0059 ± 


(0.0369 ±0.0103)?, 


V6 


= -1.2640 ± 


0.0248 


- (0.0389 ± 0.0384)?' 




v 7 = 


0.7980 ±0.0009 - 


(0.0185 ±0.0010)?, 




= 0.7880 ±0 


0128- 


(0.0268 ± 0.0202)?-, 


for / (1710) 


n" — 
v 3 — 


-1.2708 ±0.0053- 


- (0.0095 ±0.0034)?, 


v'l 


= 0.7884 ±0 


.0031 - 


(0.0041 ±0.0010)?. 



IV. CONCLUSIONS 



Before discussing the results of performed analysis, let us note that for calculating the wide resonance parameters 
(masses, total widths and coupling constants with channels) it is vital to use the poles on those sheets on which the 
poles arc not shifted due to a coupling of channels because these poles respect positions of zeros on the physical sheet. 
These appropriate sheets are numerated by 2 l (? = 1, • • ■ , N is the number of channel), i.e. II, IV, VIII,- • • . In this 
work we demonstrated this principle on the basis of analytic continuations of the S-matrix elements to unphysical 
sheets in the two-channel case. The general case of N channels can be found in other papers 0, 

It appears that neglecting the above-indicated principle can cause misunderstandings. This concerns especially the 
analyses which do not consider the structure of the Ricmann surface of the S'-matrix. For example, in literature there 
is a common opinion (delusion) that the resonance parameters should be calculated using resonance poles nearest 
to the physical region. This is right only in the one-channel case. In the multi-channel case this is not correct. It 
is obvious that, e.g., the resonance pole on sheet III, which is situated above the second threshold, is nearer to the 
physical region than the pole on sheet II from the pole cluster of the same resonance since above the KK threshold the 
physical region (an upper edge of the right-hand cut) is joined directly with sheet III. Therefore, the pole on sheet III 
influences most strongly on the energy behaviour of the amplitude and this pole will be found in the analyses, not 
taking into account the structure of the Riemann surface and the representation of resonances by the pole clusters. 
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In our model-independent approach using the uniformizing variable, we analyzed data on isoscalar S-wave processes 
tttt — > tttt, KK including the very precise NA48/2-Collaboration 7T7r-data in a threshold region. Moreover, for the tttt 
scattering the alternative data were taken: these are the data by B. Hyams et al.(1973) (set I) and by R. Kamihski 
et al.(2002) (set II) which are considerably different in the 0.76-GeV region and above 1.45 GeV. 

It was shown that when analyzing only the tttt scattering from set I, the good description of the data is achieved 
(x 2 /NDF s» 1.07) with parameters of resonances (Table HJ mainly coinciding with the ones cited as estimation of the 
PDG l| (where for the /o(600) the found pole on sheet II coincides practically with the one at around 450 — «275 MeV 
which was found in the recent dispersive 7T7r-scattering data analyses [26ll27]). Exceptions are only the mass of /o(980) 
(1001 MeV against 980±10 MeV of PDG) and the width of / (1500) (336 MeV against 109±7 MeV of PDG), though 
the most analyses, devoted to the tttt scattering, found also the state /o(980) above 1 GeV (e.g. [22j). As to the 
Jo (1500), we analyzed also the 7T7r-scattering assuming two states (narrow and wide) in this region. Description is of 
the same quality as in the first case. Parameters of the narrow Jo (1500) coincide with those preferred by the PDG. 

However, first, the satisfactory description still does not mean that it is the adequate description. The point is 
that the negative phase-shift in the background arises already on the tttt threshold. This is denoted as a pseudo- 
background. It appears to compensate for a too fast rise of the phase-shift of the amplitude, which is induced by 
the parameters of the /o(600), i.e., it indicates that these parameters are incorrect. Especially, the non zero negative 
phase shift is in contradiction with the expectation that in our parametrization the phase-shift in the tttt background 
below the KK threshold is practically zero [15| because the left-hand branch-point at s = 0, which gives a main 
contribution to the tttt background below the KK threshold, is included explicitly in the uniformizing variable. Other 
possible contributions of the left-hand cut from exchanges of the lightest mesons - the p- meson and the /o(600) - 
practically obliterate each other because vector and scalar particles contribute with the opposite signs due to gauge 
invariance. 

Second, a description of the process tttt — > KK with the resonance parameters obtained in the analysis of only 
the tttt scattering is satisfactory only for the phase shift which is due to the fact that we approximate the left-hand 
branch-point at s = 4(m|- — m 2 ) in 5 12 and £22 by the pole of the fourth power and that the pole clusters of resonances 
are chosen correctly. The module of the ^-matrix element is described satisfactory only from the KK threshold up 
to the energy about 1.15 GeV as it should be due to the two-channel unitarity (see eqs. ©)■ Above this value of 
energy the module is not described well even qualitatively. 

To this point let us also note results of our previous work Q for the coupling constants of the /q mesons with various 
channels. Despite of a preliminary character of these results, one can draw some conclusions about, e.g., the /o(600) 
and /o(980). These states turn out to have large coupling constants with the KK and especially rjrj systems, i.e., 
studying these states we deal with a multi-channel problem. Even if these states can not decay into the rjrj channel, 
their large coupling with the rjrj system should manifest itself in exchanges in the irrj scattering. 

It was shown in the two-channel approach to the tttt scattering, that the combined analysis of the coupled processes 
- the tttt scattering and tttt — > KK - is needed. This analysis was done for two sets of data (sets I and II) including 
two resonances in the 1500-MeV region. Then both above-indicated important remarks, related to the analysis only 
of the 7T7r-scattering, are ruled out. In the combined analysis the parameters of the /o(600) were changed considerably 
with the new values closer to those obtained in our previous three-channel analysis [Jj. It was shown for the data of 
sets I and II, that in the region below 1 GeV, there are two solutions, A and B, related to the cr- meson// (600) with 
the mass about 0.65 GeV and width about 0.8 GeV in the case A and m a ss m p and width about 1 GeV in the case 
B. This accords with the Weinberg prediction done on the basis of mended symmetry [30j |. Moreover, this is also in 
agreement with a refined analysis using the large- N c consistency conditions between the unitarization and resonance 
saturation suggesting m p — m a = 0(N~ 1 ) [33l |. 

Following a tradition, we speak here on the masses and total widths of resonances though the broad multi-channel 
states are represented more correctly by the pole clusters which are their mo del- independent characteristics (see the 
discussion in Sec. Ill) whereas the masses and widths are very model-dependent for wide resonances. Values of masses 
are necessary, e.g. for the mass relations of multiplets. 

The obtained values for the tttt scattering length for the A-solutions are in accordance with predictions of ChPT 
(non- linear realization of chiral symmetry), whereas the values for the B-solutions agree with predictions of chiral 
theory with linear realization of chiral symmetry (models of NJL type) . Generally, considering only description of the 
analyzed processes, it is impossible for now to prefer any of these solutions. The B-solutions for set I and II describe 
the data slightly better, whereas the obtained tttt scattering lengths for the A-solutions have more acceptable values. 
However, if one considers the problem of precise determination of the tttt scattering length a[j to be solved taking into 
account the results of the NA48/2 Collaboration [24| and the DIRAC experiment [3l[ at CERN, then the A-solutions 
should be chosen. Therefore, our final conclusion is that the agreement of our approach with ChPT and data for the 
a$ tttt scattering length favors to the A-solutions for the masses and widths of the scalar resonances. 

Finally, let us stress that our method is developed under the first principles such as analyticity, unitarity and 
Lorentz invariance, and therefore, it is free of any suppositions on dynamics except for an obvious statement that 
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a main model-independent contribution of resonances is given by the pole clusters and possible remaining small 
(model-dependent) contributions of resonances can be included in the background. 
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